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Abstract 

We consider the extended Hubbard model containing intrasite and in- 
tersite Coloumb interactions. This model we apply to a dimer as an inter- 
isting and nontrivial example of two interacting ions which possesses exact, 
analytical solution. We find the eigenvalues Ea and eigenvectors \Ec,) of 
the dimer and we represent each part Ea\Ea){Ea\ of the dimer Hamil- 
tonian (a = 1,2,..., 16) in the second quantizations with the use of the 
Hubbard and spin operators. This procedure gives a review of all compet- 
itive, intrinsic interactions together with their exactly calculated coupling 
constants expressed by the model parameters. These competitive interac- 
tions, deeply hidden in the orginal form of the extended Hubbard model, 
make an evidence that the model is extremely complex because it de- 
scribes a competition between magnetism and superconductivity. Among 
competitive interactions we can find ferromagnetic and antiferromagnetic 
interactions (also in the form appearing in the t — J model), hopping of 
the Cooper pairs between different dimer lattice sites (similar as in Kulik- 
Pedan, Penson-Kolb models), as well as, intersite Cooper pair interactions. 
We plot several coupling constants of these interactions vs Coulomb intr- 
asite interaction U to show that the competition between them strongly 
depends on the model parameters. The thermodynamical activity of each 
particular interaction, belonging to a given energy level, depends, how- 
ever, on the occupation of this level. When e.g. a given dimer energy level 
is empty the corresponding intrinsic interaction belonging to this level is 
completely passive (irrelevant). From the presented review of interactions 
it is evident that the extended Hubbard model is capable to describe the 
properties of superconductors with magnetic ordering (including also high- 
Tc superconductors) where a strong competition between magnetism and 
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superconductivity takes place. 

PACS: 71.10.-w, 71.10.Ca, 71.10.Fd, 71.20.Be, 71.27.+a, 71.28.+d 

1 Introduction 

It is commonly believed that the Hubbard model and its extentions can widely be 
used to explain many of physical phenomena observed in quite different areas of 
the solid state physics: magnetic and transport properties of transitions matals, 
their compounds and alloys, including insulator-metal transitions (cf e.g. Refs 
[1]-[13] and papers cited therein), superconductors (magnetic superconductors 
(see e.g. Refs [14], [15]), high-Tc superconductors (negative U models, cf. e.g. 
Ref. [16]; t — J based models (cf Refs [17]-[26])), fluctuating valence systems 
and heavy fermions (Anderson-like models, Ref. [27], cf also e.g. Refs [11], 
[28]), liquid He^ (cf e.g. Refs [29]-[31]), fuUcrcns (cf e.g. Refs [32]-[34]), etc.. 
To demostrate a high degree of complexity of the extended Hubbard model we 
consider one-band electronic system described by the Hamiltonian 

The indices enumerate the lattice points (Ri,Rj), is the hopping inte- 
gral, U (positive or negative) denotes the effective intrasitc Coloumb interaction, 
J^^^ and J^^^ (generally, not necessary equal) decribe the effective intersite inter- 
actions, resulting from the original Coloumb repulsion modified by the polaronic 
effects (sec e.g. Ref. [16] for details). The operators Ci^a-{cfg-) arc the electronic 
annihilation (creation) operators in the lattice site i with spin index a =t, i 
and nj,o- = cf^Ci^„. The model (1) cannot be solved exactly in a general case. 
There is, however, a special but nontrivial case of two interacting ions (a dimer 
problem) which possesses exact, analytical solution. The Hamiltonian of the 
dimer has the form (see (1)) 

Hd = -iE(ct<TC2,fT + cj^ci^<^) + U{ni,^ni^i +n2,^n2,i) 

a 

where t = — fi,2 = —^2,1- We can exactly show that Hd incorporates many 
competitive interactions, deeply hidden in its orginal form. To bring them all 
into light we will follow the method, applied earlier to the Hubbard dimer alone 
(see Ref. [35]). The application of this method to the extended Hubbard model 
for a dimer seems to be very important because it gives a review of all intrinsic 
interactions, possible to appear in this more general case. According to the 
Ref. [35] we first find the exact solution of the dimer problem (see Sect. 2). 
This solution consists of 16 energy levels Ea and corresponding eigenvectors 
\Ea) (a = 1, 2, 16). In the next step we use the equivalent expression for the 
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(1) 



(2) 



dimer Hamiltonian 



16 

Hd = ^ EaPa 

a=l 



(3) 



where Pa = \Ea){Ea\- For each a — 1, 2.., 16 the oprator Pa can be represented 
in the second quantization where we introduce Hubbard and spin operators. 
When using the relation (3) we can decompose the dimer Hamihonian (2) into 
10 parts cohecting the terms belonging to the same energy level. Such decompo- 
sition according to the exact dimer energy levels possesses two main adventages. 
First, we bring into light all competitive interactions, possible to appear within 
the exactly solvable dimer. We find ferromagnetic and antiferromagnetic inter- 
actions. Cooper pair hopping terms (similar as in Kulik-Pedan, Penson-Kolb 
models (cf. Rcfs [36], [37])) with positive and negative coupling constants, as 
well as, intcrsitc Cooper pair interactions (see e.g. Ref [16]). In this way we ex- 
actly show that the extended Hubbard model for a dimer is, in fact, a mixture of 
different intrinsic interactions leading to magnetic and superconducting proper- 
ties of the model which compete together. Second, the presence of each intrinsic 
interaction is stricly ascribed to a given energy level which can be occupied or 
empty (it depends on the model parameters and assumed average occupation 
number of electrons n G [0,4]). In other words, we can foresee which intrinsic 
interaction can be " thermodynamically" active or not (it depends on the posi- 
tion of the chemical potential ^ with respect to a given energy level) . Coupling 
constants connected with each type of intrinsic interaction are calculated exactly 
and several of them we plot as a functions of the intrasite Coulomb interaction U 
to visualize the competition between them. The method applied in this paper is 
exact in comparison with approximative approaches (perturbation expansion or 
canonical transformations) used in the literature till now (see Refs [38]- [53] and 
[11] and [12] for a review). The knowledge about possible intrinsic interactions 
within the extended Hubbard model seems to be very important because this 
model is widely used in different areas of the solid state physics. The presence of 
ferromagnetic and antiferromagnetic interactions competing with different types 
of interactions leading to superconductivity unambigously suggests the use of 
the model to magnetic superconductors (cf e.g. Refs [14], [15] and papers cited 
therein). There is, however, one important and unsolved problem which arises 
from our exact calculations for a dimer. When apply the model to a real lattice 
(real materials) we are always forced to make aproximations which can destroy 
a delicate harmony of the model leading to the overestimation or underestima- 
tion of some important competitive interactions. Thus, the problem is how to 
formulate a resonable approximation which treats all of them on equal footing. 
Besides, a high degree of complexity of the extended Hubbard model, explicitly 
shown here for a dimer, can also be a reason why the understanding of pairing 
in the high- Tc superconductors is such a difficult problem till now. It seems to 
be evident that the t — J model alone (also present in our approach), introduced 
in Refs [17]-[26] (see also Ref. [53] for a review), does not contain enough in- 
gredients to describe high-Tc superconductors because we additionally have a 
competition between ferromagnetism and antiferromagnetism accompanied by 
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a challenge between the hopping of the Cooper pairs between different lattice 
sites (similar as in Refs [36], [37]) and intersite Cooper pair interactions. 



2 Intrinsic interactions 

The eigenvalue problem for the dimer Hamiltonian (2) can be solved exactly 
when start from the vectors |ni,|, ni^j^; n2,|, n2,j.) {ui^cr = 0, 1; i = 1, 2; ct =t, J,) 
forming the Fock basis (cf. also Refs [35, [54], [55]) 



|0) = |0,0;0,0), 

|11) = |1,0;0,0), 
|12) = |0,1;0,0), 
|13) = |0,0;1,0), 
|14) = |0,0;0,1), 



|21) = 
|22) = 
|23) = 
|24) = 
|25) = 
|26) = 



1,1;0,0), 
1,0; 1,0), 
1,0;0,1), 
0,1; 1,0), 
0,1;0,1), 
0,0; 1,1), 



131) 
132) 
|33) 
|34) 



|0,1;1,1), 
|1,0;1,1), 
|1,1;0,1), 

|i,i;i,o). 



(4) 



|4) = |1,1;1,1). 



The basis vectors (4) have the form |n/3) where n = J2i a- ^i,ci= 0> 1' 2, 3, 4). The 
second index /3 (if necessary) enumerates the vectors belonging to the subspace 
of a given n. The exact solution of the dimer eigenvalue problem H£)\Ea) = 
Ea\Ea) can be obtained with the use of a standard procedure. We obtain 



El 



Eo = -t; 



= 0; l^i) 

\E2) 



111) 



|13)), 



Es=t; |i?3) = -i=(|ll)-|13)), 
EA = -t; |£;4) = ;^(|12) + |14)), 
E,=t; |i;5) = ^(|12) - |14)), 



E6 = 




1^6) = 


7l(|23) + |24)), 




Er = 


U; 


\E7) = 


;)5(|21)-|26)), 




Eg = 




\Es) = 


ai(|21) + |26))-a2(|23) 


- |24)), 


Eg = 




\E9) = 


a2(|21) + |26)) + ai(|23) 


- |24)), 


Elo = 


= J«; 


\Eio) = 


= |22), 




Ell = 


= JW; 


\Eii) = 


= |25), 





(5) 



i;i2 = t + C/ + J(l) + J(2); 

Ei3 = -t + U + J^^'^+J^^^; 

£;i4 = f+C7 + J(l)+ J(2); 

£;i5 = -t + C7 + JW + J(2); 



|i?i2) = ;^(|31) + |33)), 
|i?i3) = ^(|31)-|33)), 
|i?i4) = ^(|32) + |34)), 
1^15) = ^(|32) - |34)), 



£;i6=2(t/+ JW+J(2)). \Eie) = \4) 
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where 



C=\[^^] +4i2, (6) 



1 / (C/-J(2)) 

1 r iu - j(2)) 

Lets us introduce the Hubbard operators 

^i,<7 = Ci^aTli-a- (10) 

and spin operators 

= ^KT-na) = ^KT-<a). (11) 

•S'," = c+^c,,| = a+^ai,t (13) 

where n?^ = at^a,,,^ (i = 1, 2; cr =T, i). 

To obtain all intrinsic interactions, deeply hidden in the original form of the 
extended Hubbard model for a dimer (2) we apply the equivalent form the dimer 
Hamiltonian (3). Each product EaPa in the formula (3) where we insert Ea and 
\Ea) from the formulae (5) can be represented as a linear combination of the 
basis vectors (4) which, in turn, with the use of the Hubbard and spin operators 
(9)-(13) can be rewritten in the second quantization (see Refs. [35], [54], [55]). 
We collect all the terms together which correspond to the same energy level (we 
take into account the degeneration of the levels as e.g. E4 = E2, £"5 = Es, 
Ell = Eio, Ei4 = E12, Ei5 = E13). In this way we can split the Hamiltonian 
(2) written in the form of the formula (3) into 10 terms, corresponding to 10 
different dimer energy levels (see (5)). We obtain 

10 

Hn=J2H^'^ (14) 

i=l 

where 
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~|Z]K<^"2,<7+aJ,<.ai,a], (15) 



(T 

H^^^ = EePe = -J^'^[S^^ ■ - ^] + ^ • + Sr ■ S2+) (17) 



H^^ = ErPr=^[n\il-n^,-^)+n',{l-n1-'^)] 



1 ((72 _ fj(2)')2\ 



+ 



(6) _ .„.„_J_:^_.:^^!fc:^_2^ilr^.^_!^i 



1 (-[72 _ fj(2))2) 

+ 1^-^-- ]KT<i"2.i^2,T+&^,T4,i"ia^i.T] 

cr i=l 

h{J) = EioPio + i^iiPii = 2 J(i) • + ^] (21) 
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ff(^) - F P +p p _ (-^ + ;7 + jw + j(^)) , , 

+ti±i^±|iil±^EK>..+»i>..i. (23) 

and n"''' = J^a'^lt^ Ko- = K<^h<^ = '^i.^^i -a (« = 1,2), « = 1 when i = 2 
and i = 2 when i = 1. As a test for the correctness of the decomposition (14) we 
can sum up all the terms (15)-(24) and wc obtain the dimcr Hamiltonian in its 
original form (2). The decomposition of the dimer Hamiltonian (14) according 
to the dimer energy levels where each part is given by the expressions (15)-(24) 
has, however, very important advantages. First of all it vizualizes in an explicit 
and exact way all possible intrinsic interactions with corresponding coupling 
constants. Second, due to the decomposition according to the dimer energy 
levels we can easily see which interaction can be active or passive. It is simply 
connected with the occupation of a given level. If a dimer level is occupied 
by the electrons the intrinsic interaction ascribed to this level will be active, 
otherwise passive (it can be neglected) . The occupation, in turn, depends on the 
model parameters, temperature and the assumed average number of electrons 
— J2i ai''^i,'^) which determines the position of the chemical potential and 
decides which energy levels will be occupied or empty. 

The most important intrinsic interactions (leading to magnetism or super- 
conductivity), presented in the expressions (15)-(24) can be devided into two 
classes. First of them belongs to magnetic interactions (ferromagnetic, antiferro- 
magnetic - it depends on the sign of the parametrs J^^^, J^^^). Such interactions 
are present in the formulae (17) and (21) and describe Ising type interactions 
with coupling constants generated by J^^-* and J^^\ The formnale (17) contains 
also the transverse interaction between spins (see also later). The Heisenberg 
type magnetic interactions can be seen in the first terms of the formulae (19) 
and (20), generated by a more complex coupling constants, expressed by the 
model parameters J^^\ U and t. It is interesting to note that when J^^^ = 
the first term in the formulae (19) and (20) describes ferromagnetic or antiferro- 
magnetic interactions (it depends on the sign of U), similar to well-known t — J 
model because the coefficient ^ w ^ for large | U \ (see the formula (6) when 
= 0). The first terms in the expressions (19) and (20) can also be considered 
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Figure 1: Plot of the coupling constant Aig (the coefficient at the first term 
in the formula (19)) vs Coulomb interaction U for several values of J^^^. The 
parameter t = 1/6 for all Figs. 

as generalized t — J interactions, valid in a more general case of the extended 
Hubbard model. The coefficients Aig = + [ _ (^j^g 

phng constant at the first term in (19)) and A20 = —^3- ~ \-~ — ^^c'^ — ' ~ TT^ 
(the coupling constant at the first term in (20)) are depicted in Fig.l and Fig. 2, 
respectively, as functions of the intrasite Coulomb interaction U for serveral 
values of the parameter J^^^ . We see that Aig is negative and A20 is positive. 
Thus, the first terms in (19) and (20) describe competitive ferromagnetic and 
antiferromagnetic Heisenberg interactions. 

Let us note that the terms like 5^ |a^|a2. 162,1 = c^|C^|C2.j.C2,| (5^|aJ|ai.|5i,| = 
cJ|cJj^ci,|Ci,|) describe the hopping of the Cooper pairs (see the formulae 
(9) and (10)). Such a competitive hopping is present in the second terms 
in (18), (19) and (20). The plots of the corresponding coupling constants 

Big = l[U + C + (^^-^■^''"y") ^ (the second term in the formulae (19)) and 

B20 = j[U -C - ' ' ] (the second term in the formulae (20)) are de- 

picted in Fig.3 and Fig.4, respectively. The coupling constant Big is positive and 
i?20 negative. Thus, the second terms in the formulae (19) and (20) describe the 
hopping of the Cooper pairs with positive and negative coupling constants. Such 
terms are present in the Kulik-Pedan, Penson-Kolb superconductivity models 
(cf. Refs [36], [37]). 

The transversal products 51^5*2" (5'i~S'2^) are present in the second term 
in the formulae (17), as well as, in the scalar products of spin operators Si-S2 = 
Si''S2'' + ^iStS~+S{S}) in the first terms of the formulae (19) and (20). When 
use the second quantization representation we obtain S1S2 = — c^^c^, ci,j^C2,f 
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(5^ ^2 = — c^|Cj^|C24Ci,t) and we can see that such products of spin opera- 
tors describe the interaction of the intcrsitc Cooper pairs because we can write 

3^82 = -d2^idi^2 {Si 82^ = -^^^2^2,1) where di,2 = Ci,j.C2,t ((^2,1 = C2,j.ci,-r). 
Thus, the second term in the formulae (17) and the transversal parts in the 
first terms of the formulae (19) and (20) describe, in fact, the competitive in- 
tersite Cooper pair interactions leading to superconductivity (cf. e.g. Ref [16] 
and papers cited therein). Prom the other side the application of the resonat- 
ing valence bond approach (cf. Rofs [17]-[26]) allows to rewrite the terms like 
(5*1 • 5*2 — "^^'^ ), present in the first terms of the expressions (19) and (20), when 
introducing the following operators /2,i = ■^(c2,|Ci_| — C2,|Ci^|) and J21 = 
■^(C]^|cJl — C]^|C^|). For example, in the case of half-filling (cf. Ref. [21]) we 
can write (5i ■ S2-\) = - f2,ih,i- 

3 Conclusions 

Using the exact solution of the extended Hubbard model for a dimer we have 

reviewed the main, important interactions leading to magnetic and supercon- 
ducting properties of the model. In this way we have demonstrated that the 
relative simplicity of the extended Hubbard model is only an illusion. We have 
shown in an exact way that the magnetic properties (ferromagnetism, antifer- 
romagnetism) compete together. There is also additionally a challenge between 
two types of Cooper pair interactions (Cooper pair hopping and intersite Cooper 
pair interactions). Thus, at last, magnetic and superconducting tendencies of 
the model compete together. The final, thermodynamical properties of the ex- 
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Figure 3: Plot of the coupling constant Big (the coefficient at the second term 
in the formula (19)) vs Coulomb interaction U. The values of J''^-'(as in Figs 1 
and 2) have no significant influence on the plot in this scale. 

tended Hubbard model are then a result of this competition. Which interaction 
wins in such a rally it strongly depends on temperature (there is possible a 
change of the leader with varying temperature), model parameters and on the 
value of the assumed average number of electrons, determining the chemical 
potential of the system. Such indications can certainly be helpful when use 
the extended Hubbard model for real lattices. We have also indicated that the 
extended single-band Hubbard mock^l is, in fact, capable to describe properties 
of such superconductors where the competition between magnetism and super- 
conductivity can explain many astonishing properties of these materials (cf e.g. 
Refs [14]- [25]), including also high-Tc materials and recently synthetized C0O2 
layer based systems (cf Ref. [26] and papers cited therein). 

The nonperturbative formalism presented in this paper is, in principle, ap- 
plicable to more complicated clusters consisting e.g. of one central ion and z its 
nearest neighbours. We, however, know (see e.g. Refs [56]-[59] and papers cited 
therein) that the mathematical problems in this case exponentially increase with 
the size of the considered cluster. Another examples of a cluster approach to 
Hubbard-like models are given in Refs [60]- [65]. 
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